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(Dated: January 6, 2009) 

First measurements of the difTerential cross sections d^cr/(dp^dj/^dy''''') for the inclusive produc- 
tion of a photon in association with a heavy quark (6, c) jet arc presented, covering photon transverse 
momenta 30 < < 150 GeV, photon rapidities < 1.0, jet rapidities \y^°^\ < 0.8, and jet trans- 
verse momenta > 15 GeV. The results are based on an integrated luminosity of 1 fb~^ in pp 
collisions at ^/s = 1.96 TeV recorded with the DO detector at the Fermilab Tevatron Collider. The 
results are compared with next-to-leading order perturbative QCD predictions. 
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Photons (7) produced in association with heavy quarks 
Q (= c or b) in the final state of hadron-hadron interac- 
tions provide valuable information about the parton dis- 
tributions of the initial state hadrons 0, [1| • Such events 
are produced primarily through the QCD Compton-like 
scattering process gQ — s- jQ, which dominates up to pho- 
ton transverse momenta (pj) of ~ 90 GeV for j + c + X 
and up to ~ 120 GeV for j + b + X production, but also 
through quark-antiquark annihilation 99 — > 71? — > ^QQ- 
Consequently, "i+Q+X production is sensitive to the 6, c, 
and gluon {g) densities within the colliding hadrons, and 
can provide constraints on parton distribution functions 
(PDFs) that have substantial uncertainties 0, 01 ■ The 
heavy quark and gluon content is an important aspect of 
QCD dynamics and of the fundamental structure of the 
proton. In particular, many searches for new physics, e.g. 
for certain Higgs boson production modes [1, 0, 0, S] j will 
benefit from a more precise knowledge of the heavy quark 
and gluon content of the proton. 

This Letter presents the first measurements of the in- 
clusive differential cross sections a / [dplpdy'^ dy^°'^) for 
7 + 6 + X and ^ + c + X production in pp collisions, 
where y^' and are the photon and jet rapidities 
The results are based on an integrated luminosity of 
1.02 ± 0.06 fb-i collected with the DO detector [ul 
at the Fermilab Tevatron Collider at ^/s = 1.96 TeV. 
The highest pr (leading) photon and jet are required to 
have \y'^\ < 1.0 and < 0.8, and transverse momen- 
tum 30 < p:^, < 150 GeV and pi^' > 15 GeV. This se- 
lection allows one to probe PDFs in the range of parton- 
momentum fractions 0.01 ^ x < 0.3, and hard scatter 
scales of 9 X 10^ < = {pl^f < 2 x 10^ GeV^. Dif- 
ferential cross sections are presented for two regions of 
kinematics, defined by y'^y^°^ > and y'^y^^^ < 0. These 
two regions provide greater sensitivity to the parton x 
because they probe different sets of xi and X2 intervals, 
as discussed in Ref. [l^ . 

The triggers for this analysis identify clusters of large 
electromagnetic (EM) energy, and are based onp^ and on 
the spatial distribution of energy in the photon shower. 
The trigger efficiency is R:!96% for photon candidates with 
pI = 30 GeV and rises to nearly 100% for p^^ > 40 GeV. 



To reconstruct photon candidates, towers lj| with 
large depositions of energy are used as seeds to create 
clusters of energy in the EM calorimeter in a cone of ra- 
dius n = 0.4, where 7^ = ^(Ar;)^ -|- (A0)2 li]. Once 
an EM energy cluster is formed, the final energy -Eem 
is defined by a smaller cone of TZ = 0.2. Photon candi- 
dates are required to be isolated within the calorimeter, 
and must also have > 96% of their energy in its EM 
section. We require the sum of the total energy inside 
a cone of 7^ = 0.4, after the subtraction of £^em, to be 
< 7% of Eem- We also require the width of the energy- 
weighted shower in the most finely segmented part of 
the EM calorimeter to be consistent with that expected 
for an electromagnetic shower, and the probability for 



any track spatially matched to the photon EM cluster 
to be <0.1%. Background from dijet events containing 
TT*^ and 7] mesons that can mimic photon signatures is 
also rejected using an artificial neural network for iden- 
tifying photons (7- ANN), described in Ref. 1^. The 
requirement that the 7-ANN output be > 0.7, combined 
with all other photon selection critera, reduces the di- 
jet event efficiency to 0.1-0.5%. We calculate photon 
detection efficiencies using a Monte Carlo (MC) simula- 
tion. Signal events are generated using PYTHIA [lij and 
processed through a GEANT-based [15| simulation of the 
detector geometry and response, and reconstructed using 
the same software as for the data. The MC efficiencies are 
calibrated to those in data using small correction factors 
measured in Z — > e^e~ samples. The total efficiency of 
the above photon selection criteria is 63-80%, depending 
on The systematic uncertainties on these values are 
5%, and are mainly due to uncertainties in the isolation, 
the track-match veto, and the 7-ANN requirements. 

At least one jet must be present in each event. Jets 
are reconstructed using the DO Run II algorithm [l^ 
with a radius of 0.5. The efficiency for a jet to be re- 
constructed and to satisfy the jet identification criteria 
is 93%, 96.5%, and 94.5% for light (m, d, s quark or g), 
c, and b jets at p^ = 30 GeV and increases to « 98% at 
p^ = 150 GeV, independent of the jet flavor. The impact 
from uncertainties on jet energy scale, jet energy resolu- 
tion, and difference in energy response between light and 
6(c) jets is found to be between 8%(6%) and 2%(2%) 
for p!^^ between 15 GeV and 150 GeV. The leading jet is 
also required to have at least two associated tracks with 
Pt > 0.5 GeV and the track leading in px must have 
Pt > 1-0 GeV, and each track must have at least one hit 
in the silicon microstrip tracker. These criteria ensure 
that the jet has sufficient information to be classified as 
a heavy-fiavor (HF) candidate. Light jets are suppressed 
using a dedicated artificial neural network (6- ANN) [l3| 
that exploits the longer lifetimes of heavy-fiavor hadrons 
relative to their lighter counterparts. The leading jet is 
required to have a 6-ANN output > 0.85. Depending on 
p"^, this selection is 55-62% efficient for 7 -I- 6 jet, and 
11-12% efficient for 7 -I- c jet events, with 3-5% relative 
uncertainties on these values. Only 0.2-1% of light jets 
are misidentified as heavy-fiavor jets. 

A primary collision vertex with >3 tracks is required 
within 35 cm of the center of the detector along the beam 
axis. The missing transverse momentum in the event is 
required to be < 0.7p^ so as to suppress background from 
cosmic-ray muons and W ^ decays. Such a require- 
ment is highly efficient for signal, achieving an efficiency 
> 96% even for events with semi-leptonic heavy-flavor 
quark decays. 

About 13,000 events remain in the data sample after 
applying all selection criteria. Background for photons, 
stemming mainly from dijet events in which one jet is 
misidentified as a photon, is still present in this sample. 
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FIG. 1: Distribution of observed events for PHF-jot after all 
selection criteria for the bin 50 < < 70 GeV. The distribu- 
tions for the b, c, and light jet templates are shown normalized 
to their fitted fraction. Error bars on the templates represent 
combined uncertainties from statistics of the MC and the fit- 
ted jet fiavor fractions, while the data contain just statistical 
uncertainties. Fits in the other bins are of similar quality. 



FIG. 2: The j + b + X and 'y + c+X differential cross sections 
as a function of in the two regions y~'y''^^ > and y^'y^^^ < 
0. The uncertainties on the data points include statistical 
and systematic contributions added in quadrature. The NLO 
pQCD predictions using CTEQ6.6M PDFs are indicated by the 
dotted lines. 



To estimate the photon purity, a template fitting tech- 
nique is employed [3]. The 7- ANN distribution in data 
is fitted to a linear combination of templates for pho- 
tons and jets obtained from simulated 7 + jet and dijet 
samples, respectively. An independent fit is performed in 
each p7p bin, yielding photon purities between 51% and 
93% for 30 < < 150 GeV. The fractional contribu- 
tions of h and c jets are determined by fitting templates 
of PHF-jct = -Inll^^track the data, where P^^^^ is 
the probability that a track originates from the primary 
vertex, based on the significance of the track's distance of 
closest approach to the primary vertex. All tracks within 
the jet cone are used in the fit, except the one with low- 
est value of Ptrack- Jets from b quarks usually have large 
values of -PuF-jct, whereas light jets mostly have small 
values, as their tracks originate from the primary ver- 
tex. Templates are used for the shape information of the 
^HF-jct distributions. For h and c jets these are extracted 
from MC events whereas the light jet template is taken 
from a data sample enriched in light jets, which is cor- 
rected for contributions from b and c quarks. The result 
of a maximum likelihood fit, normalized to the number of 
events in data, is shown in Fig. [1] for 50 < < 70 GeV. 
The estimated fractions of 6 and c jets in all p7p bins vary 
between 25-34% and 40-48%, respectively. The corre- 
sponding uncertainties range between 7-24%, dominated 
at higher p7p by the limited data statistics. 

The differential cross sections are extracted in five bins 
of p'^ and in the two regions of y'^y-'®*, and are all listed 
in Table m The measured cross sections are corrected for 
the effect of finite calorimeter energy resolution affecting 



pip using the unfolding procedure described in Ref. [20t . 
Such corrections are 1-3%. The measured differential 
cross sections are shown in Fig. [5] for 7 + & + X and 
^ -\- c + X production as a function of pTp for the jet 
and photon rapidity intervals in question. The cross sec- 
tions fall by more than three orders of magnitude in the 
range 30 < < 150 GeV. The statistical uncertainty 
on the results ranges from 2% in the first pTp bin to « 9% 
in the last bin, while the total systematic uncertainty 
varies between 15% and 28%. The main uncertainty at 
low pTp is due to the photon purity (10.5%) and the heavy- 
flavor fraction fit (9%). At higher pTp^ the uncertainty is 
dominated by the heavy-flavor fraction. Other significant 
uncertainties result from the jet-selection efficiency (be- 
tween 8% and 2%), the photon selection efficiency (5%), 
and the luminosity (6.1%) [l3|. Systematic uncertainties 
have a 60-68% correlation between adjacent pTp bins for 
30 < < 50 GeV and 20-30% for pT^ >70 GeV. 

Next-to-leading order (NLO) perturbative QGD 
(pQGD) predictions, with the renormalization scale ^ij, 
factorization scale /i p , and fragmentation scale pf, all set 
to pTp, are also given in Table U and compared to data in 
Fig. O These predictions [19] are are based on techniques 
used to calculate the cross section analytically [IJ], and 
the ratios of the measured to the predicted cross sections 
are shown in Fig. [3) 

The uncertainty from the choice of the scale is es- 
timated through a simultaneous variation of all three 
scales by a factor of two, i.e., to fJ-R,Fj = 0.5bZ and 
2pJ. The predictions utilize CTEQ6.6M PDFs [J], and 
are corrected for effects of parton-to-hadron fragmenta- 
tion. This correction for b (c) jets varies from 7.5% (3%) 
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TABLE I: The "/ + b + X and j + c + X cross sections in bins of in the two regions y^'y''^^ > and y''y''^^ < together with 
statistical, Sastat, and systematic, Sasyst, uncertainties. The theory cross sections athcory are taken from Ref. p^ . 



y'-'y'"*' > ^^yJ"' < 





bin 


(Pt) 


Cross section i 


^stat 




^ theory 


(Pt) 


Cross section Sastat 




•^theory 




(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 


(pb/GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 


(pb/GeV) 


J + b + X 


30-40 


34.1 


2.73x10'^ 


1.5 


18.5 


2.96xl0~' 


34.1 


2.23x10"' 


1.6 


19.1 


2.45x10"' 




40-50 


44.3 


1.09x10"^ 


2.5 


15.5 


9.31x10"^ 


44.2 


9.53x10"^ 


2.6 


16.0 


8.18x10"^ 




50-70 


57.6 


2.72x10"^ 


3.3 


15.2 


2.66x10"^ 


57.4 


2.67x10"^ 


3.3 


15.3 


2.22x10"^ 




70-90 


78.7 


6.21x10"^ 


6.6 


20.8 


6.39x10"^ 


78.3 


6.10x10"^ 


6.7 


20.8 


5.49x10"^ 




90-150 


108.3 


1.23x10'^ 


8.2 


26.2 


1.11x10"^ 


110.0 


1.09x10"^ 


8.9 


25.7 


1.05x10"^ 


7 + c + X 


30-40 


34.1 


1.90 


1.5 


18.1 


2.02 


34.1 


1.56 


1.6 


18.7 


1.59 




40-50 


44.3 


5.14x10"^ 


2.5 


17.7 


5.82x10"^ 


44.2 


4.51x10"' 


2.6 


18.1 


4.56x10"' 




50-70 


57.6 


1.53x10"^ 


3.3 


17.9 


1.41x10"^ 


57.4 


1.50x10"' 


3.3 


18.0 


1.10x10"' 




70-90 


78.7 


4.45x10"^ 


6.6 


21.3 


2.85x10"^ 


78.3 


4.39x10"^ 


6.7 


21.3 


2.22x10"^ 




90-150 


108.3 


9.63x10"^ 


8.2 


27.5 


3.69x10"^ 


110.0 


8.57x10"^ 


8.9 


27.0 


3.28x10"^ 



h4 

f1.6| 

S1.4- 
n 

°1.2: 
1: 
0.8: 
0.6; 
0.4 ; 
0.2^ 

3.5 T 
3 

2.5 T 
2 ': 

1.5 
1 T 

0.5 



D0, L , = 1 .0 fb 

' ml 

y > 
y+b + X 



1/1 < 0.8 
|y|<1-0 
p!^'>15GeV 



Y + b + X 




data / theory 

CTEQ6.6M PDF uncertainty 
ICBHPS/CTEQ6.6M 
IC sea-like /CTEQ6.6M 
Scale uncertainty 



yy"' > 
y + c + X 



yY"' < 
y + c + X 



40 60 80 100 120 140 40 60 



80 100 120 140 

p;(GeV) 



FIG. 3: The data-to-theory ratio of cross sections as a func- 
tion of for 7 + b + X and 7 + c + X in the regions 
y7yjet y Q y''y'°^ < 0. The uncertainties on the data in- 
clude both statistical (inner line) and full uncertainties (entire 
error bar). Also shown are the uncertainties on the theoretical 
pQCD scales and the CTEQ6.6M PDFs. The scale uncertain- 
ties are shown as dotted lines and the PDF uncertainties by 
the shaded regions. The ratio of the standard CTEQ6.6M pre- 
diction to two models of intrinsic charm is also shown. 



at 30 < < 40 GeV to 1% at 90 < pj. < 150 GeV. 

The pQCD prediction agrees with the measured cross 
sections for ^ + b+X production over the entire range, 
and with -f+c+X production for < 70 GeV. For > 
70 GeV, the measured ^ + c + X cross section is higher 
than the prediction by about 1.6-2.2 standard deviations 
(including only the experimental uncertainties) with the 
difference increasing with growing . 

Parameterizations for two models containing intrinsic 



charm (IC) have been included in CTEQ6.6 0, and their 
ratios to the standard CTEQ predictions are also shown 
in Fig.[3l Both non-perturbative models predict a higher 
^ + c+X cross section. In the case of the BHPS model Q 
it grows with pj. The observed difference may also be 
caused by an underestimated contribution from the g — > 
QQ splitting in the annihilation process that dominates 
forpj, > 90 GeV 



In conclusion, we have performed the first measure- 
ment of the differential cross section of inclusive pho- 
ton production in association with heavy flavor (b and 
c) jets at a pp collider. The results cover the range 
30 < < 150 GeV, |?/'^| < 1.0, and \y^°^\ < 0.8. The 
measured cross sections provide information about b, c, 
and gluon PDFs for 0.01 < a; < 0.3. NLO pQCD predic- 
tions using CTEQ6.6M PDFs jl§| for 7 -I- 6 -I- X production 
agree with the measurements over the entire pj. range. 
We observe disagreement between theory and data for 
J + c + X production for pj, > 70 GeV. 
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